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An outstanding issue in the treatment of boundaries in general relativity is the lack of a local 
geometric interpretation of the necessary boundary data. For the Cauchy problem, the initial data 
is supplied by the 3-metric and extrinsic curvature of the initial Cauchy hypersurface.. This Cauchy 
data determines a solution to Einstein's equations which is unique up to a diffeomorphism. Here, 
we show how three pieces of boundary data, which are associated locally with the geometry of the 
boundary, likewise determine a solution of the initial-boundary value problem which is unique up to 
a diffeomorphism. One piece of this data, constructed from the extrinsic curvature of the boundary, 
determines the dynamical evolution of the boundary. The other two pieces constitute a conformal 
class of rank-2, positive definite metrics, which represent the two gravitational degrees of freedom. 
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I. INTRODUCTION 

There exists a well-posed Cauchy problem for Einstein's equation which has the important property of geometric 
^ uniqueness, i.e. local geometric data representing the 3-metric and extrinsic curvature of the initial Cauchy hypersur- 
(2J[) face determines a spacetime metric metric gab which is unique up to diffeomorphism. There are two formulations of 
'—'I the initial-boundary valueproblem (IBVP) which are strongly well-posed, the Friedrich-Nagy formulation [T] and the 
^ . harmonic formulation 0, Q, but neither give rise to a such a local version of geometric uniqueness. For the harmonic 
^ formulation, there exists a nonlocal (in time) version in which the geometric interpretation of the boundary data 
, depends upon a background metric constructed from the initial Cauchy data As a result, the boundary data have 
' a geometric interpretation which is nonlocal in time. In this work, we show how boundary data for the gravitational 
00 , field can be posed which are locally determined by the geometry of the boundary, in the same sense as the Cauchy 
C5 • data. 

' In a Cauchy problem, initial data on a spacelike hypersurface determine a solution in the domain of dependence 
C3 ■ 2-'(5o) (which consists of those points whose past directed characteristics all intersect iSq). In the IBVP, data on a 
' timelike boundary T transverse to Sq are used to further extend the solution to the domain of dependence I?(5o U T). 
. Strong well-posedness @ guarantees the existence of a unique solution which depends continuously on both the initial 
^ ' data and the boundary data. 

The primary application of the gravitational IBVP is the simulation of an isolated astrophysical system containing 
neutron stars and black holes. The standard approach in numerical relativity, as in computational studies of other 
hyperbolic systems, is to introduce an artificial outer boundary T, which is coincident with the boundary of the 
" ■ computational grid and whose cross-sections are spheres surrounding the system. The ability to compute the details 
of the gravitational radiation produced by compact astrophysical sources, such as coalescing black holes, is of major 
importance to the success of gravitational wave astronomy. If the simulation of such systems is not based upon a 
strongly well-posed IBVP then the results cannot be trusted in the domain of dependence of the outer boundary. For 
comprehensive reviews of the gravitational IBVP see d, 0| • 

Geometric uniqueness has the practical application of allowing numerical simulations with the same initial data 
but carried out with different formulations and different gauge conditions to produce geometrically equivalent (dif- 
feomorphic) spacetimes. For a recent discussion of this issue, see [8]. For hyperbolic systems which are stable under 
lower order perturbations, the global solution in the spacetime manifold Ai can be obtained by patching together 
local solutions, i.e. the problem can be localized. Thus, for purposes of treating the underlying geometrical nature of 
the boundary data, it suffices to concentrate on the local problem in the neighborhood of a point on the boundary. 
That is the approach taken in this paper. 

In the Friedrich-Nagy formulation, there are three essential pieces of boundary data that have geometrical or phys- 
ical significance. One is the trace K of the extrinsic curvature Kab of the boundary, which geometrically determines 
the location of the boundary. (Note that the coordinate specification of the location of the boundary is pure gauge 
information and does not determine its location geometrically in the same sense that a curve is determined geomet- 
rically by its acceleration, given its initial position and velocity.) Two other pieces of data in the Friedrich-Nagy 
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formulation, which are related to the gravitational radiation degrees of freedom, are supplied by a combination of the 
Weyl tensor components and ^"4, in the Newman-Penrose notation Q. The remaining boundary data specify the 
gauge freedom. 

The Friedrich-Nagy formulation is based upon a symmetric hyperbolic Einstcin-Bianchi system, with evolution vari- 
ables consisting of an orthonormal tetrad, the associated connection coefficients and the Weyl curvature components. 
Although it differs from the metric based formulations used in numerical relativity, the requirement of three pieces of 
geometric boundary data should be universally applicable. (Statements found in the literature that only two pieces of 
boundary data suffice to specify the physical or geometrical properties of the gravitational field are misleading. They 
are only true when the boundary has been geometrically specified, e.g. for an r = const boundary in a background 
Schwarzschild geometry.) 

The outgoing null vector K"" and ingoing null vector used in defining and 5*4, respectively, are determined 
by the unit normal N"" to the boundary and a choice of unit timelike vector tangent to the boundary according to 

/ST'^ = r° + iV" , L^^T^-N"". (1.1) 

The choice of T° represents gauge freedom in this data. Friedrich and Nagy are careful to point out that this gauge 
freedom prevents interpreting ^fp and ^'4 as purely geometric data. 

This shortcoming could perhaps be avoided by choosing these null vectors to be principle null directions of the 
Weyl tensor. However, in a general spacetime this would lead to four possible choices which would then have to be 
incorporated (in some yet unkown way) into a well-posed problem. An alternative, proposed in Q , is to base the data 
on the eigenvectors V"" determined by the trace-free part of the extrinsic curvature according to 

where Hab is the intrinsic 3-metric of the boundary. For a spherical worldtube in Minkowski space, this picks out 
a locally preferred timelike direction T". This suggests that the approach might extend to a suitably round outer 
boundary of an isolated system. However, it is again not clear whether such an approach can be properly incorporated 
into the evolution system. 

Here we consider geometric boundary data for metric based formulations of the IB VP. Our main result is that, 
along with the initial Cauchy data, a geometrically unique spacetime is determined by three pieces of boundary data 
related locally to the intrinsic metric Hab and extrinsic curvature Kab of the boundary. More specifically, the boundary 
data consist of a conformal class of rank-2, positive definite metrics, which represent the two gravitational degrees of 
freedom, and an associated component of the extrinsic curvature, which determines the dynamical evolution of the 
boundary. The rank-2 metric on the 3-dimensional boundary picks out a timelike eigen-direction, with null eigenvalue, 
which, up to normalization, provides a geometric choice of unit timelike vector tangent to the boundary. In Sec.llll 
we discuss the underlying geometry and present our main result as a Local Geometric Data Theorem. 

The demonstration that this data lead to a strongly well-posed IBVP is carried out using the harmonic reduction 
of Einstein's equations to ten wave equations, as was the method used in establishing the analogous result for the 
Cauchy problem fl^. In doing so, the three pieces of local geometric boundary data must be supplemented by seven 
additional boundary conditions. Four of these conditions are supplied by the harmonic coordinate conditions. The 
other three fix the remaining freedom in the choice of harmonic coordinates. In Sec. IIIIl the resulting harmonic 
IBVP is reduced to a set of partial differential equations in the frozen coefficient formalism, which are subject to a 
combination of Dirichlet and Neumann boundary conditions. 

The demonstration that the strong well-posedness of the frozen coefficient version of the harmonic IBVP extends 
to the full quasilinear problem was given in 0, Q for the case of Sommerfeld boundary conditions. In Sec. llVi we 
demonstrate how this result can be extended to Dirichlet and Neumann conditions. 

The strong well-posedness of the frozen coefficient harmonic IBVP with local geometric boundary data is then 
treated in Sec.|Vl The key idea is that the full set of boundary conditions can be applied sequentially, similar to the 
approach followed in 0, Q except now applied to a set of Dirichlet and Neumann conditions rather than Sommerfeld 
conditions. Our main result is then established in Sec. IV Dl 

When the emphasis is on geometric issues we use abstract tensor indices, e.g. v"" to denote a vector field, and 
when the specific spacetime coordinates — {t,x^) are introduced we use the corresponding coordinate indices, e.g. 
v'^ = (w*, w'). 

II. THE INITIAL-BOUNDARY DATA 

We begin with a review of the initial data for the Cauchy problem. The standard treatment of the Cauchy problem 
introduces a time foliation St, with future directed unit normal n^. The embedding of St in the spacetime manifold 
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A4 then gives rise to the decomposition of the spacetime metric 

gab = -naUb + hab, (2.1) 

where hab is the 3-mctric intrinsic to St ■ Geometric initial data are determined by the intrinsic metric hat and extrinsic 
curvature kab = hf^S/ c^b of the initial Cauchy hypersurface iSq, where Vq is the covariant derivative associated with 
gab- These data are subject to the Hamiltonian and momentum constraints 

2G''^nanb =11- kabk"^ + fc^ = 0, (2.2) 
hlC^Ua^ Dbik"'' - h'^^'k) = 0, (2.3) 

where TZ is the curvature scalar and Db is the covariant derivative associated with hab- 

The remaining initial data necessary to determine a unique spacetime metric consist of gauge information, i.e. data 
that affect the solution only by a diffeomorphism. In the 3 + 1 formulation of Einstein's equations, the gauge freedom 
in the metric is governed by the choice of an evolution field 

r = an" + /3^ p^ua = 0, (2.4) 

with lapse a and shift /S". The lapse relates the unit future-directed normal to the time foliation, according to 

na = -aVat- (2.5) 

The evolution field is transverse but not in general normal to the Cauchy hypersurfaces so that it determines the shift 
according to 

= h^tK (2.6) 

The initial data required for the formulation of a well-posed Cauchy problem depends upon the choice of hyperbolic 
reduction of Einstein's equations. Here we consider the hyperbolic reduction associated with harmonic coordinates, 
as used in the classic work of Choquet-Bruhat Generalized harmonic coordinates = {t,x^) = {t,x,y,z) are 
functionally independent solutions of the curved space scalar wave equation 

g'^'VaybX^' = ~f (2.7) 

where t^^{g,x) are harmonic gauge source functions [ll| . Thus the harmonic coordinates are determined by initial 
data x'' = (0, a;') on So, subject to the initial conditions dixf^ = S^. 

In terms of the connection coefficients T^^^, the harmonic coordinate conditions are 

:= - f ^ = 0, (2.8) 

where 

1 



^ gP^r^ ^ ^ d^^gP^l g^detg,,. (2.9) 
The hyperbolic reduction of the Einstein tensor results from setting 

E'"' := - W^T'> + ^^g^^VpCP = 0, (2.10) 

where is treated formally as a vector field in constructing the "covariant" derivatives V^C". 

When the harmonic conditions (|2.8p are satisfied, the principle part of (|2.10l) reduces to a curved space wave 
operator acting on the densitizcd metric, i.e. 

Et"" = -^=g°''^dadB{y^gf"') + lower order terms = 0. (2.11) 

Thus the harmonic evolution equations (j2.10p are quasilinear wave equations for the components of the densitized 
metric \/~gg^'^ - The well-posedness of the Cauchy problem for the harmonic system p.lOp follows from known results 
for systems of quasilinear wave equations [13]. Such results are local in time since there is no general theory for the 
global existence of solutions to nonlinear equations. 
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Constraint preservation results from applying the contracted Bianchi identity V i^^G^'^ = to (|2.10l) . which leads to 
the homogeneous wave equation 



V'VpC^ + i?^C'' = 0. (2.12) 
If the initial data enforce 

Cl5o=0 (2.13) 

and 

dtC'Uo = (2.14) 

then C = is the unique solution of (|2.12[) . It is straightforward to satisfy (I2.13|) by algebraically determining the 
initial values of dtg^^ in terms of the initial values of g^^^ and their spatial derivatives. In order to see how to satisfy 
(I2.14P note that the reduced equations (I2.10[) imply 

G^^^n^ = n^V^^C") - ^ri^VpC. (2.15) 

As a result, if 

G^^n,\s,=Q, (2.16) 

i.e. if the Hamiltonian and momentum constraints are satisfied by the initial data, and if the reduced equations (j2.10p 
are satisfied, then 

KV^^r^-in^VpCll^o =0. (2.17) 

It is easy to check that if C^l^^ =0 then (^17)) implies (pHl) . 

By standard results, the Hamiltonian and momentum constraints on the initial data, along with the reduced 
evolution equations (|2.10l) . imply that the initial conditions (|2.13p and (I2.14p required for preserving the harmonic 
conditions are satisfied. Conversely, if the Hamiltonian and momentum constraints are satisfied initially then (j2.15p 
ensures that they will be preserved under harmonic evolution. Thus the conditions = substitute for the constraints 
of the generalized harmonic formulation. 

These results show that the free initial gauge data in the harmonic formulation consist of the initial values of 
the lapse and shift. For simplicity, we set the initial lapse to unity and the initial shift to zero so that the metric 
components satisfy 

5" = -l, .9" = 0, t = Q. (2.18) 

Along with the initial geometric data /i*-' and A;*^ , these determine a unique solution to the Cauchy problem in the 
harmonic gauge. In geometric terms, hab and kab determine a solution which is unique up to a diffeomorphism. 

We now formulate the additional geometric data necessary for the IBVP. In the IBVP, there is another natural 
decomposition of the metric at the boundary T, 

gab = NaNb + Hab, (2.19) 

where Na is the unit outward normal and Hab is the 3-metric intrinsic to T ■ The boundary T intersects the initial 
Cauchy hypersurface at an edge Bq. In general, the spacelikc normal Na to T is not orthogonal to the timelike 
normal Ua to Sq. The geometric initial data must now also include the hyperbolic angle Qq at the edge given by 

sinheo = iVa"lBo- (2-20) 

The initial velocity of the boundary with respect to the initial Cauchy hypersurface is determined by Qq. 

We represent the local geometric data on the boundary T which encodes the two gravitational degrees of freedom 
by a conformal class {Qab} of rank 2 positive definite, symmetric tensor fields, i.e. the class defined by the equiv- 
alence relation Qab = ^^Qab, > 0. Here {Qab\ is identified with the intrinsic metric of the boundary by further 
decomposing (12.19^ into 

Hab — —TaTb + Qab, (2-21) 
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where Qab belongs to {Qab} and QabT^ — 0. Here T° is the future directed unit timelike vector tangent to the 
boundary which is geometrically picked out as an eigen-vector of {Qab\ with with null eigenvalue. On the initial 
Cauchy hypersurface, we identify Qab with the intrinsic metric of the edge Bq. 

The remaining part of the local geometric data on the boundary T is obtained from its extrinsic curvature 

Kab = HIV.N^. (2.22) 

In the Friedrich-Nagy formulation of the IBVP, the trace K = W^^Kab forms part of the boundary data. Using the 
fact that Hab has ( — t- +) signature, Friedrich and Nagy show that when K is expressed in terms of a boundary 
defining function it gives rise to a wave equation for that function which geometrically determines the location of the 
boundary. For the methods we use here to establish the strong well-posedness of a metric formulation of the IBVP, 
there does not appear to be a way to incorporate K into the boundary data. However, the alternative component 

L = {H"^ - T''T^)Kab, (2.23) 

which is geometrically determined by Qab, does supply the data in the necessary form. Because {H°'^ — j,arpb^ ^^^^ 
has ( — h +) signature, L geometrically determines the location of the boundary by the same construction used by 
Friedrich and Nagy. 

We can know state our main result. 

Local Geometric Data Theorem:. Cauchy data hab cmd kab on Sq along with edge data Qq on Bq and boundary 
data {Qab} and L on T determine a solution of the vacuum Einstein equations (locally in time) which is unique up 
to a diffeomorphism. All data are assumed to be smooth and compatible. 

Here the Cauchy data must satisfy the Hamiltonian and momentum constraints but the boundary data are constraint 
free subject to compatibility with the Cauchy data, e.g. the restriction of {Qab} and hab to Bq must lead to conformally 
equivalent 2-metrics. Together {Qab} and L supply three pieces of local geometric boundary data. 

As for the case of the Cauchy problem, additional data, which control the gauge degrees of freedom, are necessary 
to determine a unique solution. This gauge data depend upon the particular hyperbolic reduction used to formulate 
the IBVP. In the formulation of a strongly well-posed harmonic IBVP, the Einstein equations reduce to 10 wave 
equations for the components of the metric, so that 10 boundary conditions are necessary. In addition to the 3 pieces 
of geometric data {Qab} and L, the harmonic conditions (|2.8I) supply 4 boundary conditions, as described in Sec. IIIII 
Thus 3 more pieces of gauge data on the boundary are necessary to specify completely the harmonic coordinate 
freedom. This data pins down the values of the harmonic coordinates on the boundary. 

A non-zero value of the hyperbolic angle 9o presents a technical complication in prescribing these 3 pieces of 
harmonic gauge data. However, the value of &o can be adjusted to zero by carrying out a Cauchy evolution in the 
neighborhood of Bq to a new choice of 5o, which keeps Bq unchanged. Since the Cauchy problem is well-posed, the 
initial data for this modified problem depend continuously on the initial data for the original problem. Consequently, 
the original IBVP is strongly well-posed if the IBVP for the modified problem with 0o = is strongly well-posed. In 
the following, we assume that this has been carried out. (Otherwise, the technical details in constructing a convenient 
gauge for establishing a well-posed IBVP become more complicated; cf. jSj] where the case 6o 7^ is treated.) The 
requirement that Oq = implies 

A^a^lfjo = SO that r"|B„ = (2.24) 

Since harmonic coordinates are solutions of the curved space scalar wave equation, they are determined by the initial 
data and boundary data for a scalar wave. The boundary data for these coordinates can be specified in any form 
which leads to a strongly well-posed IBVP. For our present purpose, we consider Dirichlet or Neumann boundary data. 
In order to investigate the possible choices, let = (t, x, x^) be Gaussian normal coordinates in the neighborhood of 
the boundary, where x^ = {ij, z), with the boundary at 2; = for i > 0. In these coordinates, the metric has the form 

gf.odx^dx'' = dx^ + Hj jdx^dx'^ x' = (f, x^). (2.25) 

In addition, we choose T°'daX^ = on the boundary so that = 0, i.e. T" has vanishing shift relative to the 
i-foliation of the boundary. Then the intrinsic 3-metric of the boundary has components 

Hjjdx'dx-^ = TfT^dP + Q^^dx^dx^ . (2.26) 

(Note that the lapse intrinsic to the boundary could also be set to 1 locally, but globally in time such Gaussian 
coordinates for the intrinsic boundary metric could introduce focusing singularities.) 
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Harmonic coordinates x^^ = {t,x,x"^), = {y,z), can now be introduced by solving an IBVP for the scalar wave 
equation (|2.7p . On the boundary we prescribe the Dirichlet data x = x = 0, so that the boundary is given by x = 0. 
For the remaining harmonic coordinates we prescribe the Neumann data 

dx"^ dt 

— =0, 7TT=0, a; = 0, 

ox ox 

so that on the boundary 



^ dx dx" ^ dx dx^ 

Similarly g^^ — on the boundary, which is consistent with the initial condition p.24p at the edge B{). In summary, 
the boundary freedom in the choice of harmonic coordinates allows us to set 

5"*lr = 5"^lr = 0, x|r = 0. (2.27) 

III. REDUCTION TO PDES 

In order to reduce the IBVP to a set of partial differential equations (PDEs) with the initial-boundary data described 
in Sec. |TT1 we express the harmonic Einstein equations (12. lip in the form 

g^Pd^dpi^ggn^F^". (3.1) 

where the forcing F^'^ represents lower order terms which do not enter the principal part. Since the harmonic gauge 
source functions play no essential role, we set = 0. 

In the harmonic coordinates defined in Sec.|lll the initial data at t = 0, with the gauge conditions p.l8|) . consist of 

g^l = 3*^=0, .9" = -l, 

dtg'^ ^ -h'^, dtiV^g'') ^ -d,{Vhh'^), dt{V^g'')^0. (3.2) 
The boundary data at x = consist of the geometric data 

Q'^^ L = {Q-'' ~2T''T'')Kab^-^V9^{Q'"' -2T''T'')dM (3.3) 

for some representative Q°'^ = il~^Q"'^ of the conformal class {Q"*"}, along with the boundary gauge data (I2.27p . 

3"'=0, g^^^O. (3.4) 

In addition we enforce the harmonic constraint on the boundary, 

d^.{V^gn\x=o^O. (3.5) 

We now formulate the PDEs for the frozen coefficient version of the problem. The material in See's. IIVI and Sec. IVl 
shows that the strong well-posedness of this frozen coefficient problem extends to the quasilinear problem. 

In this approach, the problem is localized in the neighborhood of a point p on the boundary and the wave operator 
in (j3.1l) is frozen to its value at p, 

g'^^lx.d^dp. (3.6) 

By a constant linear transformation of the harmonic coordinates which keeps the z-direction fixed, we can then set 
5"^ Up = v"^ (ths Minkowski metric). In doing so, the x-direction remains aligned with N"^ and we can further align 
the ^-direction with T°. Then in the neighborhood of p we linearize the equations about the Minkowski metric. In 
terms of these coordinates x'^ = {t, x, x"^) = {t, x, y, z) and the linearized variable 

= ^93^'' (3.7) 

the system ()3.ip takes the frozen coeffcient form 



dx dx^ - - dx dx 
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= F, x>0, t>0, (3.8) 
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with forcing matrix F. 

In these coordinates, Q*^ and the conformal metric data consist of 

qAB ^ q1/2qAB^ q ^ (3 9) 

In the hnearized approximation, this reduces to 

and the extrinsic curvature data reduces to 

L = iQ^'" - 2T''T'')Kab - -^a,(27^^ +72'^' + 7^^). (3.11) 

Here it is assumed that this boundary data is sufficiently close to that for a boundary with vanishing extrinsic 
curvature in a flat spacetime to allow the iterative construction of a solution to the quasilinear problem. However, 
this can be arranged by considering the conformally transformed metric g'^^^ — A^g^jy, where A << 1 is a positive 
constant; cf. p. 262 of [12]. Then L' = XL and, in the stretched coordinates x'^ — xt^ + \^^{x^ ~Xp)i the transformed 
metric has components 17^/^/(2;') = g^v{x) = "qf^v + 0(A). 

The boundary conditions for this linearized system are 



i(7?/?^-7-) = gi(t,y,z), 



7-* = 0, 



7^2^ = 0, 
7- = 0, 

5*7*" + 5.7"" + dyj^y + 9.7"" = 0, 
5t7*^ + 5.7"^ + dyj^y + 9,7^' = 0, 
5t7*' + 9.7"" + dyjy^- + d.Y'- = 0, 
9*7" + 9.7*" + dy^'y + 9,7*^ = 0. 



3.12) 

3.13) 

3.14) 

3.15) 
3.16) 
3.17) 
3.18) 
3.19) 
3.20) 
3.21) 



The Dirichlet data qi and q2 are determined from the two conformally invariant degrees of freedom Q^^ . The 
Neumann data qs are determined from the extrinsic curvature scalar component L. The Dirichlet conditions p.lSp 
- p.l7p arise from the boundary conditions (13. 4p on the harmonic coordinates. The boundary conditions p.lSp - 
p.2ip arise from the harmonic constraints p.Sp . 



IV. ENERGY ESTIMATES FOR QUASILINEAR WAVE PROBLEMS WITH SOMMERFELD, 
DIRICHLET OR NEUMANN BOUNDARY CONDITIONS 

We establish the strong well-posedness of quasilinear IBVP for wave equations with Dirichlet and Neumann bound- 
ary conditions by an approach similar to that carried out in [3] for Sommerfeld boundary conditions. We begin by 
reviewing the Sommerfeld case. 
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A. Sommerfeld boundary conditions 

The energy estimates in Sections 1-4 of 3 established that the the frozen coefficient version of the harmonic IBVP 
with Sommerfeld boundary conditions is unique and depends continuously on the data. In Appendix A of Q, this 
result was extended to the strong well-posedness of the quasilinear problem. Here we first consider a slightly simplified 
version of the problem treated in Q. We show that local existence theorems and energy estimates for second order 
quasilinear wave equations can be obtained in the same way as for first order symmetric hyperbolic systems. It 
all depends on a priori estimates for arbitrarily high derivatives of the solutions of linear equations with variable 
coefficients. 

Consider the half-plane problem 

Uu = Pu + Ru + F, X > 0, < y < oo, (4.1) 

with Sommerfeld- type boundary conditions at a; = 0, 

a{ut + 7w) — Ux + q , a>0, 7>0 strictly positive constants, (4.2) 
smooth boundary data q{t, y) and smooth initial data 

uiQ.x.y) = fi{x,y), ut{0,x,y) ^ f2{x,y). (4.3) 
Here the subscripts (i, x, y, z) denote partial derivatives, e.g Ut — dtu, 

Pu = {aux)x + {buy)y - 2^ut - 

and 

Ru = CiUt + C2Ux + C^Uy + C4U. 

Ru are terms of lower (first and zeroth) differential order. Also, we use the notation 

{u,v), = (u,^); {u,v)b, \\u\\%^ {u,u)b 

to denote the L2 scalar product and norm over the half-plane and boundary, respectively. 

All coefficients and data are smooth real functions and a > oq > 0, 6 > &o > 0, where ao,6o a-i'e strictly positive 
constants. The initial data are compatible with the boundary conditions. Here 7 > is a constant obtained by the 
change of variables u — ?> e^^u' and then deleting the "prime". This introduces the term 7^||u|p in the energy 

E \\utf + {ux, aux) + {uyMy) + 7^l|wf , (4.4) 

which provides an estimate of 

Lemma. There is an energy estimate which is stable against lower order perturbations. 
Proof: Integration by parts gives 

dtE = dt^WutW"^ + (ux.aux) + {uy.buy) +-f'^\\u\\'^) 

= -47||ut|p +2(ut,F) +2(ut,i?M) - 2{ut,aUx)B 
+ at\\ux\\^ + bt\\uy\\^ 

< const{E + \\F\\^)~2{ut,aUx)B- (4.5) 
Using the boundary conditions, we obtain 

-{ut,aUx)B = -{ut,aaut)B - {ut,aaju)B + {ut,aq)B < -^aQajdt\\u\\% + const(||u||| + \\q\\%)- 
Therefore (|4.5p implies 

dt{E + aoa7||u|||) < const(ii; -I- \\Ff + \\u\\l + Ml)- (4.6) 

This proves the lemma. 
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Now we can estimate the derivatives. Let v = Uy, w = Uf. Differentiation of the differential equation gives 

Vtt = PV + RV + RyU + {ayUx)x + {byV)y + Fy, 

Wtt = Pw + Rw + Rtu + {atUj.)x + (htv)y + Ft. (4.7) 

Here RyU and Rtu are hnear combinations of first derivatives of u which we have aheady estimated and can be 
considered part of the forcing. 

The differential equation (14. ip tells us that 

O'Uxx — wt ~ bvy + terms we have already estimated. 

Thus is lower order with respect to v and w and, except for lower order terms, v and w are solutions of the 
same differential equation as u. They obey the same boundary conditions with data y) and qt(t, y), respectively. 
Therefore we can estimate all second derivatives. Repeating the process, we can estimate any number of derivatives. 

We can now proceed in the same way as in Q, where we have considered first order systems to obtain existence 
theorems for equations with variable coefficients. We approximate the differential equation by a stable difference 
approximation and prove, using summation by parts, that the corresponding estimates for the divided differences 
hold independently of gridsize. In the limit of vanishing gridsize, we obtain the existence theorem. Since we can 
estimate any number of derivatives, it is well known, using Sobolev's theorem, that we can obtain similar, although 
local in time, estimates for quasilinear systems. By the same iterative methods as for first order symmetric hyperbolic 
systems it follows that strong well-posedness extends locally in time to the quasilinear case, as well as other standard 
results such as the principle of finite speed of propagation. 

Remark. There are no difficulties to extend the results to three spatial dimensions. 



B. Homogeneous Dirichlet and Neumann conditions 



If we replace the Sommerfeld boundary conditions by homogeneous Dirichlet or Neumann conditions, then the 
estimates in Sec. IIV Al hold with boundary data q = 0. 

Now we consider the half-plane problem for wave equations with inhomogeneous Dirichlet or Neumann boundary 
conditions. We want to show that we can transform these problems into problems with homogeneous boundary 
conditions by changing the forcing and the initial data. As a model problem, we consider the half-plane problem 

utt = (a(t,x,y)ux)^ + (b{t,x,y)uy)y + F{t,x,y), (4.8) 
a; > 0, — oo < y < oo, t > 0, 

with initial conditions 

u{0,x,y) = fi{x,y), Uf {0,x,y) ^ f2{x,y), (4.9) 

and Dirichlet boundary conditions 

uit,0,y)^q{t,y). (4.10) 
We assume that all coefficients and data are compatible and smooth. We make a change of variable 

u{t, X, y) = uit, X, y) - (p{x)q{t, y). (4.11) 
Here Lp{x) is a smooth function, with (^(0) = 1, which decays exponentially. Then 

{t(t, 0, y) = 0, i.e. -u satisfies homogeneous Dirichlet boundary conditions. (4-12) 

By dm]), 

Utt = Utt - i(p{x)q(t,y))tt, 
{aux)x = {aux)x - iai(pix)q{t,y)x)^ , (4-13) 

{bUy)y = {bUy)y~ {b{(p{x)q{t,y)y)y. 
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Finally, by (j4.8|) . (|4.12l) and (|4.13p we obtain the differential equation with modified forcing term 

iitt = {a{t, X, y)u^)^ + {b{t, x, y)uy)y +F + F, (4.14) 

which satisfies homogenous Dirichlet boundary conditions. By assumption, is a smooth function and F is composed 
of a, 6, Lp and q and their first two derivatives. Since derivatives are smooth functions, F is also a smooth function. 
Therefore u{t, x, y) satisfies the estimates arrived at in Sec. IIV Al 
Now we consider ()4.8p with the Neumann boundary condition 



u^{t,0,y)^q{t,y). (4.15) 

We make again the transformation (|4.1ip but now with tpxiO) = 1, and obtain the corresponding result. 

As an illustration of how the estimates extend to higher derivatives, consider the half-plane problem ()4.8p with a 
homogeneous Dirichlet boundary condition for a = b = 1 (which poses no restriction), 

utt ^ Uxx + Uyy + F, u{t,0,y) = 0. (4-16) 

By differentiation, we obtain 

Uytt = Uyxx + Uyyy + Fy, Uy {t , , y) = , (4-17) 

uttt = utxx + utyy + Ft, ut{t,0,y) = 0. (4-18) 

As in (14. 7p . we introduce the variables 

V = Uy, w — ut- (4-19) 

Then (|iTfl) . become 

Vtt = Vxx+Vyy + Fy, V {t , , y) = , (4.20) 

wtt = Wxx + Wyy + Ft, w{t,0,y)=0. (4-21) 
Integration by parts then gives us an energy estimate for 

WvtW + + WvyW + + Ik.f + Ik.f • (4.22) 

By (|4J6| and (|4T9)) we obtain 

UXX + F ^ Utt - Uyy ^ Wt - Vy. 

Therefore, by (|4.22l) . we obtain a bound for 

We obtain a bound for ||ita;2;a;|p in the same way by replacing v and w by 

v^^^ = Uyy, = utf (4.23) 

Now we obtain the differential equations 

vil^ = vi'J + + Fyy , {t, 0,y) = 0, 

= + wi'J + Ftt, w^'^ {t, 0, y) = 0, (4.24) 

and we obtain energy estimates for 

\\4'^r + iw^^r + H'^r + w^'^r + + , (4.25) 

which we can express in terms of u according to 

htyyf + IWxyyf + \\Uyyy\\^ + \\utttf + \\Uxtt\\^ + \\Uyttf ■ (4.26) 

By differentiation of (|4.16p with respect to x, 

'^xxx — ^xtt ^xyy Fx- (4.27) 
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From (|4.26p . we already have estimates for HitxttlP and Therefore we also obtain an estimate for Hm^xxIP- 

This process can be continued. 

Our result is not restricted to the model problem but is valid in general. For example, we can replace (|4.8|) by the 
corresponding half-plane problem in three spatial dimensions. 

Remark. In many problems, surface waves, glancing waves and other waves specific to the boundary are important. 
In that case, there is no energy estimate and the above technique does not activate these phenomena. Instead, in such 
cases, we split the problem into two problems; one with homogeneous boundary conditions and another where only the 
boundary conditions do not vanish, i.e. the forcing and the initial values are zero. The first is covered by Sec. \IV\ 
The second we treat by Fourier- Laplace techniques. For examples, see Jldi . 



V. THE STRONG WELL-POSEDNESS OF THE IBVP FOR THE HARMONIC EINSTEIN EQUATIONS 

Here we establish the strong well-posedness of the gravitational IBVP for the system (|3.8p with boundary conditions 
(|3.12p - p.2ip determined by local geometric boundary data and harmonic coordinated conditions. In order to 
illustrate how the estimates in Sec. IIVI apply we progress through a sequence of model problems. 



A. Model problem I: The harmonic Einstein equations in one spatial dimension 

First consider the half-plane problem in the frozen coefficient formalism of the harmonic Einstein equations for the 
system of wave equations in one space variable 

/ tt tx \ 

{-d^ + dl) ^..j^F, x>0,t>0, (5.1) 

with forcing matrix F. In standard notation, we treat the system in the sequential order 

(1) 9,V" = 9^*- + Fi, 

(2) 9,V" = 92^-- + F2, (5.2) 

(3) dh''^dir + Fs. 

Here 7*^(i, x), 7^^(t, x), 7**(t, x) denote the dependent variables which we want to determine on the half-plane. The 
forcing terms Fi{t,x), F2{t,x), F3{t,x) are smooth functions of (t, x). 

The solution of our problem is determined by the initial data corresponding to (|3.2p . 

the Dirichlet boundary condition 7*^(i, 0) = q{t) 

or the Neumann boundary condition d^j*^ {t,0) = q{t), (5-3) 

and the harmonic constraints, which we apply on the boundary in the sequential order 

9t7*"(i,0)-f a,7""(i,0) = 0, (5.4) 
9t7"(t,0)-)-9,7*"(i,0) = 0. (5.5) 

We start with the wave equation for 7*^ with smooth boundary data (j5.3p and smooth compatible initial data. 
By means of the transformation (|4.1ip in Sec. IIVBI we modify the forcing so that the variables satisfy homogeneous 
boundary conditions, which we denote by 

q{t) = 0. (5.6) 

Then there is an energy estimate. The problem is strongly well-posed and we can solve the wave equation to estimate 
7*^ and its derivatives on the boundary, as well as in the interior x > 0, in terms of the data. Next we use the 
constraint (|5.4[) and obtain Neumann boundary data for 7^^ from 9(7*^ (i, 0). We again use the transformation (|4.1ip 
so that dx^^^{t,0) = 0, using the notation ()5.6p . The resulting wave problem for 7^^ with homogeneous Neumann 
data is strongly well-posed so that we can estimate 7^^(t,a;) and its derivatives. Finally, we obtain the same result 
for 7**, using the constraint (|5.5p and the transformation (|4.1ip . 
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B. Model problem II: The harmonic Einstein equations in two spatial dimensions 

Now consider the half-plane problem in frozen coefficient formalism in two spatial dimensions, 

^tx ^xx ^xv\p^ X > 0, i > 0, - cx)< y < CX3, 

^ty ^xy ^yy J 

where F again represents the forcing. The components 7^^, 7*^ and 7^^ satisfy Dirichlet or Neumann boundary 
conditions. The initial data correspond to (|3.2p . 

The harmonic constraints are applied on the boundary in the sequential order 

9t7*"(i,0,y)+5,7""(i,0,y)+a,7"nt,0,2/) = 0, (5.7) 
dti'^{t,Q,v)+dxl''\t,0,y)+dy^yy{t,Q,y) - 0, (5.8) 
aa"(i,0,2/) + a.7*"(i,0,y)+a,7*ni,0,y) = 0. (5.9) 

We proceed essentially in the same way as for model problem I. We use the transformation (|4.11l) such that wave 
equations for 7**^, 7*^ and 7^^ satisfy homogeneous Dirichlet or Neumann boundary conditions. Then there is an 
energy estimate for these variables and we can use the constraints to obtain estimates for the remaining variables. 
The constraint (|5.7p determines Neumann boundary data for 7^^(i, 0,y). After using the transformation (|4.1ip . it 
reduces to 

9,7--(i,0,y) = (5.10) 

and the resulting wave problem for 7^^ is strongly well-posed. Thus we can estimate 7^^(i,x) and its derivatives. 
Similarly, the constraint (15. 8p determines Dirichlet boundary data for 7*^(t, 0, y). After the transformation (|4.1ip . the 
constraint (j5.8p reduces to 

9*7*^(^,0, y) = (5.11) 

and the resulting wave problem for 7*^ is strongly well-posed. The constraint ()5.8p now determines Dirichlet boundary 
data for 7**(i, 0, y) and we can use the transformation (|4.1ip to obtain a a strongly well-posed problem for 7**. 



C. Model problem III: The harmonic Einstein equations in three spatial dimensions 



We now consider the half-plane problem for the linearized harmonic equations in three spatial dimensions 

/7* 



^d^+dl + dl+dl 





^tx 


jty 


7 


tx 




ryXy 


i' 


ty 


ryXy 


ryyy 




tz 




^y 


1 







7 

rjty ^xy ^yy ^yz 



= F, (5.12) 



V7* 

a;>0, i>0, — 00 < y < 00, — 00 < z < cjo, 

where F represents the forcing, 7^^, 7^^, 7^^, 7*^, 7^^ and 7^^ satisfy Dirichlet or Neumann boundary conditions 
and the initial data correspond to p.2[) . 

The harmonic constraints are applied on the boundary in the sequential order 

5t7*" + 5.7"" + dyj^y + 5,7"' = 0, (5.13) 

da'^ + d^^'^y + dy-fyy + d^jy^ = o, (5.14) 

5*7*' + 5,7^" + a^7^^ + a,7^^ = 0, (5.15) 

5t7" + a.7*" + 5j,7*« + a.7*" = 0. (5.16) 

We proceed in the same way as in two space dimensions. We use the transformation (j4.11l) so that the six wave 
equations for 7^^, 7^^, 7^^, 7*^, 7^^ and 7^^ satisfy homogeneous Dirichlet or Neumann boundary conditions. Then 
there is an energy estimate for these variables and their derivatives. We then use the constraints to obtain estimates 
for the remaining variables. 

The constraints (|5.13p - (|5.15p determine Neumann boundary data for dxj^'^it, 0, y, z) and Dirichlet boundary data 
for j*y{t, 0, y, z) and 7*^(i, 0, y, z) in terms of the previously estimated variables.. After using the transformation (|4.1ip . 
the resulting wave problem is strongly well-posed so that we can estimate 7^^, 7*^ and 7*^ and their derivatives. The 
constraint (j5.16p then provides Dirichlet data 7**(t, 0, y, z) for the remaining variable in terms of previously estimated 
variables After the transformation (|4.1ip . the resulting wave problem for 7** is strongly well-posed. 
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D. The harmonic Einstein equations with local geometric data 

Now we turn to the 3-dimensional harmonic Einstein system (|5.12l) with boundary conditions p.l2p - p.2ip 
determined by local geometric boundary data and harmonic coordinate conditions, as prescribed in Sec. IIIIl After 
applying the transformation (|4.1ip . the conformal metric data (13.121) - p.l3p reduce to the homogeneous Dirichlet 
form 

(^yy_^-)(t,0,y,z) = 0, jy^-{t,0,y,z) = 0, (5.17) 

and the extrinsic curvature data L (|3.14p reduce to the homogeneous Neumann form 

d^jyy + + 27--)(i, 0, y, z) = 0. (5.18) 

The boundary gauge data (|3.15p - p.l7p are already in the homogeneous Dirichlet form 

7*"(t, 0, y, z) = ^-y{t, 0, y, z) = ^^^t, 0, y, z) = 0. (5.19) 

The remaining boundary conditions are supplied by the harmonic constraints ()5.13p - (|5.16p . 

The situation is similar to model problem III but simpler since the gauge conditions (|5.19p are already homogeneous 
and imply that the constraint (|5.13p reduces to the homogeneous form 

a,7""(i,0,y,z) = (5.20) 

so that (|5.18p reduces to 

d^{^yy +Y-'-){t,Q,y,z) = Q. (5.21) 

Together the homogeneous boundary conditions (I5.17p . (I5.19p . (I5.20p and (|5.21l) determine strongly well-posed wave 
problems for (7^^ — 7^^), 7^^, 7*^^, 7^^, 7^^^, and (7^^-1-7^^), respectively. Thus we can estimate those variables and 
their derivatives. Now we can proceed as in Model problem III to use the constraints (|5.14p - (|5.16p in sequential order 
to determine the required estimates for the remaining three variables 7*^, 7*^ and 7**. Along with the applicability to 
the quasilnear problem outlined in Sec. llVi this establishes the Local Geometric Data Theorem proposed in Sec. [Til 

VI. DISCUSSION 

We have shown that the three pieces of locally geometric boundary data supplied by a conformal class of rank-2 
metrics {Qab} and an associated extrinsic curvature component L determine, up to gauge, the necessary boundary 
data for a solution of Einstein's equations. This resolves an outstanding issue in the geometrical understanding of the 
gravitational IB VP. 

The method used to establish this result also broadens the possible formulations of a strongly well-posed harmonic 
IBVP. The method in 2, 3] based upon Sommerfeld conditions has been extended to include Dirichlet and Neumann 
conditions, subject to the sequential structure necessary to enforce the harmonic constraints. For computational 
applications, Sommerfeld conditions are most benevolent because they allow numerical error to leave the grid. It 
is therefore somewhat discordant with numerical application that a treatment of the boundary based upon locally 
geometric data must apparently include at least two Dirichlet conditions, associated with {Qab\, and one Neumann 
condition, associated with L. 

There are many options in formulating a suitable combination of Dirichlet, Neumann and Sommerfeld conditions 
for a strongly weel-posed problem, provided the sequential structure is maintained. However, none of these options 
seem to be consistent with a choice of locally geometric boundary data different from {Qab} and L. it is possible that 
other choices of might be established by a different approach. For example, had we used the trace K of the extrinsic 
curvature of the boundary instead of the component L then (j5.2ip would have been replaced by 

5.(7^^+7""-7")(i,0,2/,z) = 0. (6.1) 

This does not fit into the sequential structure for applying the constraints but it is still possible that the resulting 
system is strongly well-posed. It depends upon whether the higher derivatives can be estimated by some generalization 
of the scalar treatment given in Sec. IIVI 

An additional issue of practical importance is the formulation of a strongly well-posed IBVP for the 3-1-1 approach 
which has historically played a major role in numerical relativity [isj . In the 3-1-1 formalism, instead of the 10 wave 
equations of the harmonic system, Einstein's equations are reduced to a pair of 6 first order in time equations for hab 
and fcati, supplemented by 4 conditions which determine the lapse and shift. Perhaps the geometric insight provided 
by our results can be transferred to shed light on this problem. 
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